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Abstract 

We present a simple proof of the stability of the hydrogen molecule 

{M + M + m~ m~) . It does not rely on the proton-to-electron mass ratio M/m 

being very large, and actually holds for arbitrary values of M/m. Some 

asymmetric molecules of the type (mfm^m^m^) are also stable. Possible 

applications to molecules containing antiparticles and to exotic hadrons in the 

quark model are briefly outlined. 
13.35.Cs 
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I. INTRODUCTION 



The study of the hydrogen molecule is usually carried out in the framework of the Born- 
Oppenheimer approximation . The effective potential between the two protons exhibits a 
deep pocket of attraction, and there is no doubt that stability survives quantum fluctuations 
and other corrections due to proton mass being finite. A simple and rigourous derivation 
of the stability is, however, desirable, as well as a systematic study of other 4-unit-charge 
systems. 

A proof of the stability of the hydrogen molecule was briefly outlined in a recent com- 
munication 0. It does not rely on the proton mass M being very large as compared to the 
electron mass m. It is actually based on a result obtained for the M = m case: Hylleraas 
and Ore || have shown in 1947 that the positronium molecule (e + e + e~e~) is stable against 
dissociation into two positronium atoms. The stability of the hydrogen molecule, and more 
generally of any (M + M + m~m~) configuration, can be proved using a rescaled version of 
the wave function of the positronium molecule. 

If one aims at mathematical rigour, one can object that in the analysis of Refs. @,f|, 
it is implicitly assumed that the lowest threshold consists of two neutral atoms. It is also 
desirable to make connection with the usual approach to hydrogen binding, based on the 
Born-Oppenheimer limit M/m — > oo. In a recent letter the question of the thresholds is 
seriously addressed, and an alternative variational wave function was proposed, which proves 
stability for M > 0.144m, and reduces to the standard Heitler-London wave function in the 
limit M/m — > oo. A more detailed account of this work on the thresholds and on the 
stability for M ^> m will be published elsewhere ||. 

In the present paper, we explain how stability in the equal-mass case (e + e + e~e~) implies 
stability for all hydrogen-like configurations (M + M + m~m~). Moreover, from the calculated 
binding energy of (e + e + e~e~) with respect to its threshold (e + e~) + (e + e~), we predict a 
minimal extension of the stability region toward even more asymmetric configurations of the 
type (mim^m^m^). 
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Perhaps the stability of the Coulombic molecule (pE + f2~S~) with baryon number B = 4 
and strangeness S = 4 will be thought as being of marginal interest. (For the zoology of 
charged leptons, mesons, and baryons, we refer to the Review of Particles Properties ||.) 
There is however a very general and fundamental concern about understanding why matter 
sometimes shows up in large compounds and sometimes splits into small clusters. We hope 
to better penetrate the mechanism of molecular binding by studying how it is sensitive to 
the masses of the constituents which are involved. Moreover, our new approach to hydrogen 
binding allows for some generalisations. As briefly outlined in Sec. [VIII|, exotic molecules 
involving a mixture of matter and antimatter might well be stable, as long as annihilation 
is neglected. The stability problem also exists in other fields. In quark physics, only the 
minimal configurations have been seen so far, namely mesons made out of a quark and an 
antiquark, and baryons consisting of three quarks. We shall mention in Sec. [IX] that the 
patterns observed in molecular physics might provide some guidance in quark physics, or at 
least in simple quark models, to guess which flavour configurations are the most likely to 
host stable multiquarks. 

II. THE STABILITY PROBLEM 

Let us consider four elementary particles, with masses rrii, and unit electric charges 
qi = (+1, +1, —1, —1). There are obvious symmetries: an overall charge conjugaison, as well 
as (1 <-> 2) and (3 <-> 4) exchanges, leave the problem of stability invariant. So one can 
assume without loss of generality that 

rrii > ™>2, ™>3 > "7,4, (1) 

and that the positive charges are on the average heavier than the negative ones, if there is 
any overall asymmetry. 

In most cases, instability manifests itself by dissociation into two neutral atoms. For 
some values of the mass ratios, the lowest threshold consists of a 3-body ion and a charge, 
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for instance (777177727773) + + mj: we know that (mfm^m^) is stable in a wide range of values 
of 777.2 and 7773/7771 in particular, every ion with mi = m 2 is stable ||; moreover, if 

mi, 7772, and 7773 are large, the binding energy of (777177727773) 4 " might well exceed (in absolute 
value) that of every other threshold. However, a molecule (mfm^m^m^) would hardly 
decay into (777177727773 ) + + mj , since the long-range Coulomb attraction would keep the ion 
and the charge together. In other words, if the lowest threshold is made of a 3-body ion 
and a charge, the stability of the 4-body molecule is guarantied. We shall come back to this 
point in Sect. [\T|. 

The difficult task is thus to compare the 4-body ground state to the thresholds consisting 
of two neutral atoms. Out of the two possible arrangements, (1,3) + (2,4) and (1,4) + 
(2,3), the former is more stable, if one assumes the mass ordering ([[]). This can be seen 
directly from the Bohr formula for the 2-body energy, E 2 (mf ,m~) = — 1/(20^), with = 
m^ 1 + mj 1 being the inverse reduced mass. This is in fact a rather general result [f[0|| . If 
v (r) is an universal potential (independent of the masses), the reduced 2-body Hamiltonian 
h = ctp 2 /2 + v (r) depends on the inverse reduced mass a linearly. Then the ground state 
energy E 2 (a) is a concave function of a [fCT]] , and if ([!]) is satisfied, i? 2 (ai3) + -^2(^24) < 
£2(^14) +£ 2 («23)- 

Some states lying above the threshold are metastable, since their decay involves painful 
rearrangements and tunnelings. We shall not consider them. 

One should also mention unnatural parity states. Imagine a molecule carrying orbital 
angular momentum and parity J p = 1 + , besides the intrinsic spins and parities of its 
constituents. It cannot decay into two ground-state atoms j p = + separated by an orbital 
momentum i = 1. Its actual threshold is made of a ground state (IS) and an orbital 
excitation (2P). We shall not consider this problem further, but simply mention that some 
of the general methods and results presented in this article can be extended to unnatural- 
parity configurations. 
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III. THE EQUAL-MASS CASE 



The stability of the positronium molecule was shown in 1947 by Hylleraas and Ore || , 
who used an elegant variational method. They first got rid of the scale by noticing that if 
^(Fj) is a trial wave function, with norm and expectation values of the kinetic and potential 
energies written as 



n =(*|*), t = <tf|T|tf), v=(V\V\V), 



(2) 



then the best rescaling of the type <fi — \&(rj/A) yields a minimum 

E - 



(3) 



which corresponds to (<f)\T\(j)) = —(<f>\V\<f))/2, i.e. the same sharing of the kinetic and poten- 
tial energies as for the exact solution. This extension of the virial theorem to variational 
approximations is well known [0 . 

The frozen-scale wave function of Hylleraas and Ore contains a single parameter: 



* = exp -^(n3 + ri4 + r 23 + r 24 ) cosh ~(t-i 3 - r u - r 2 3 + r 2 4) 

Explicit integration leads to (a misprint in Ref. || is corrected below) 
33 33 - 22p 2 + 5/3 4 

T) = -I- 

16 16(1 - /? 2 ) 3 
21 3/P 21 - 6/3 2 + /? 4 
*~1F~~2~ + 8(1 - /3 2 ) 3 
19 21 - 18(3 2 + 5/3 4 



(4) 



(5a) 
(5b) 



V 6 + 4(1 - /3 2 ) 3 



1 



(1-/32)2 



5/3 2 



1 + 1 



(5c) 



4/3 4 ' 8(3 2 4/3 6 l-/3 2 _ 
to be inserted in (|J) and (H), leading to a minimum E = —0.5042 near j3 2 = 0.48. 

It is rather easy to generalise the calculation of Ref. || to a trial wave function of the 
type 



* = °i ex P - a i( r i3 + ru + r 2 3 + r 24 ) cosh bi(r 13 - r u - r 23 + r 24 ) 



(6) 



but one does not gain much [13] . As analysed for instance by Ho []14 and by the authors he 



quotes, some explicit r 12 and r 34 dependence is needed in the wave function to improve the 
accuracy. The last variational calculation gives an energy one can express as 



£(e + e + e"e~) = -(1 + x 1 )/2, x x = 0.03196, (7) 

where X\ represents the fraction of additional binding with respect to the threshold, which 
is E th = —1/2 in our units. 

IV. FROM THE POSITRONIUM TO THE HYDROGEN MOLECULE 

Once the stability of (e + e + e~e~) is established, it is extremely simple to derive that of 
every (M + M + m~m~) configuration []2|,||. Let us, indeed, rewrite the Hamiltonian as 

H = H S + H A (8a) 

H * = {-m + i) (pi +pi +Pl +Pl) + v (8 b ) 

where H$ is even under charge conjugaison, and if a odd. One notices that H and ifs have 
the same dissociation threshold, namely — (M _1 + m -1 ) -1 in our units. Now, H$ is stable 
by simple rescaling of the positronium-molecule case, and the ground state of H is lower 
than that of H§. This latter result comes from the variational principle, with the ground 
state of Hs, ^(Hs), used as a trial wave function for H: 

E(H) < (*(H S )\H S + Ha\*(H s )) = E(H S ) (9) 

since (^(ifs^-^A^^s)) vanishes, by symmetry considerations. 

It is amazing that every exotic molecule with the same structure as hydrogen is stable. 
An example is D + D + Q~Q~ , with charm C = 2 and strangeness S = —6. 

Previously, Abdel-Raouf |T(| and Rebane [F7 stressed the regularity of the binding energy 



as a function of the mass ratio M/m, but missed the fact that stability in the positronium 
case implies stability for other (M + M + m~ m~ ) configurations. 
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If one fixes for instance the scale by m = 1, one can study the lowest energy E of 
(M + M + m~m~) as a function of s = 1/M, or equivalently the fraction x(s) defined as the 
positronium case as 

l + x(s) 



i:(.\r m 



m m 



l + s 



(10) 



E th = — 1/(1 + s) being the threshold energy. The fraction x(s) takes the values X\ = 0.0303 
in the positronium case (s = 1), and xq = 0.1745 for s = 0, the limit of hydrogen with 
infinitely massive protons. The values are taken from the compilation by Rebane |17| . This 
means the inequality @, E(H) < E{H$), is actually observed. 



Since s enters the Hamiltonian linearly, E(s) is a concave function of s fLl |. On can also 
combine convexity and scaling and deduce that —1/E(s) is also concave, and this provides 



a stronger constraint [11|. The concavity of —1/E leads to an upper bound on x(s) for 
intermediate configurations 

-l 



1 + x(s) < 



1 



+ 



l + x 1 2(1 + x ) 



(11) 



V. ASYMMETRIC MOLECULES 

Let us go by steps toward (mim^m^m^), which denotes the most general case. Consider 
first a different arrangement of only two masses, (M + m + M~m~). We saw that it is stable 
for M = m. In case of a large asymmetry M ^> m, it seems to become unstable: the 
hydrogen-antihydrogen system, for instance, hardly survives its decay into (pp) + (e + e~). 

A minimal domain of stability around M = m can be derived from the variational 
principle. One can fix the scale by setting M _1 = 1 — y and mT x = 1 + y. Then 

H(M + m + M-m-) = ^ (p? + f 2 + f 3 + + V + | (-p? + f 2 - f 3 + . (12) 

The reasoning is the same as in Sec. |TV|. The last term, antisymmetric under simultaneous 
(1 <-> 2) and (3 <-> 4) exchanges, lowers the ground-state energy. Then from Eq. (|7|), 
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E(y)<E th (0)(l + Xl ). (13) 
Meanwhile the threshold becomes 

EM = -— ^— - — ^— = (14) 
{U) 4(1-2/) 4(l + y) 1-y 2 1 ; 

Thus stability remains at least as long as y 2 < 1 — (1 + xi) -1 , i.e. for 

M 

0.70 < — < 1.43 . (15) 
m 

Now, if one starts from (M + m + M~m~), and introduce four different masses rrii such that 

mi > m 3 — m <±-, m \ 1 + ^jT 1 = 2M _1 , m^ 1 + 1 = 2m~ 1 , (16) 

then one can rewrite the Hamiltonian as 

H(mim2m^m^) = H(M + m + M~m~) 

+ - m," 1 )^ - vl) + ^K" 1 - m^)(pl - pi), (17) 

and, again, the two last terms cannot do anything but lower the ground-state energy, while 
the threshold does not change. Hence stability is improved. 

For instance (fi~£ + H~p) is stable, as long as strong interactions are neglected, since the 
constituent masses are here 1.67, 1.20, 1.32 and 0.94GeV/c 2 , which fulfil flTJ) and flTJ|). We 
are not too surprised to learn that the positronium hydride (pe + e~e~) is stable, by a small 
margin |18[]: it corresponds to M/m = 2, not too far outside our minimal range flT5|) , and 
benefits from the large asymmetry between particles 1 and 3. 

The role of symmetries should be emphasised. Consider for instance the stability of 
(1 + , 1 + , l~,m~), whose threshold is 

£ th (m)=£th(l)^^ (18) 

where 5 = mT l — 1 measures the departure from the symmetric case of the positronium 
molecule. One can first split the Hamiltonian into 



S 



H{m) = H(l) + -pi 



(19) 



and say that the ground state should lie below the first order estimate. The expectation 
value of P4 is easily estimated from the virial theorem and the symmetries of the unperturbed 
wave function. Then 



E(m) < E(l) 1 - 



(20) 



leading to a sufficient condition for stability 

1 + 5/4 



1 52 



< 1 + Xt. 



(21) 



(l-5/4)(l + 5/2) 

If one instead splits the Hamiltonian into pieces of well-identified behaviour under permu- 
tations 



H = 



" / -2 1 ^2 — *2 1 ->2\ 1 " / — *2 1 ^2\ 

-Pi + P 2 - Ps + P4J + 7 (-P2 + P4J , 



(22) 



one gets 



£(m) < E(l) 



corresponding to a wider range of stability 



1 + 5/ 4 ' 



(23) 



'1+5/4 2 1 5 2 

■ r~J~^ — = H 1 < l+Xl. 

1 + 5/2 16 



(24) 



VI. VERY ASYMMETRIC MOLECULES 

We already mentioned that some configurations are qualitatively different from the 
positronium molecule, since their lowest threshold does not consist of two neutral atoms. 
These configurations cannot studied by starting from (e + e + e~e~), and implementing more 
and more asymmetry. 

Consider for instance the molecule (A+A+f2~e~), with charmed hyperons, or any similar 
case where three particles are much heavier than the fourth one, and form a stable ion. Its 
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stability results from the (A+A+f2~) ion lying below its dissociation threshold into an atom 
and a charge, as is the case for every mfm^m^ ion with m\ = m 2 |||7|]. The electron is 
bound around this ion, but the electronic energy is negligible. 

We expect interesting properties in the transition region where a threshold consisting of 
3-body ion and a charge becomes equal to the lowest threshold made of two neutral atoms. 
The competing thresholds are likely to generate some attraction in the whole 4-body system, 
and make it stable. This particular situation deserves some specific investigations. 

VII. TENTATIVE GRAPHICAL SUMMARY 

Thanks to the scaling properties of the Coulomb interaction, the study of a system 
of N given charges requires only (N — 1) independent variables to scan all possible mass 
distributions. 

In the case of 3 unit charges qi = ±(—1,1,1), one could introduce 2 independent mass 
ratios. It was however found more convenient to use 3 barycentric coordinates a^, with 
< on < 1, and a normalisation J2 a i — 1- Each case corresponds to a point inside an 
equilateral triangle, in which a stability frontier separates the stability from the instability 
areas. One can choose the masses as barycentric coordinates, as in Fig. 0. The inverse 
masses of Fig. ^| prove more suited for the mathematical analysis of the observed convex 
behaviour of the frontier ||. 

The generalisation to AT = 4 charges corresponds the inner volume of a regular tetrahe- 
dron, so that the stability frontier keeps the symmetries of the problem. We shall summarise 
our results and our guesses both with the masses as variables (J2 m i — 1), or the inverse 
masses a« oc m^ 1 [J2 a i — 1); restricting ourselves to very schematic drawings of the un- 
folded tetrahedrons. Future theoretical and numerical works will hopefully make it possible 
to determine the frontier more accurately. 

Consider first the representation in terms of the masses, in a regular tetrahedron (ABCD) 
of unit height, so that Yl m i — 1> with mi being the distance to the face (BCD), etc. This 
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corresponds to Fig. |3|. 

The summit A stands for the masses m 8 = (1,0,0,0), i.e. the limit of a positronium 
hybride (pe + e~e~) with an infinitely massive proton. It is stable, at least for rri2 = = rriA 
||18|| . This means the stability domain reaches A, at least along the symmetry axis. A 
similar situation is of course observed near B, C, and D. To our knowledge, there is no result 
available on how stability survives when the light masses rri2, and are different, in 
particular for mi ^> 772,3 ~ 1714 ^> m2, corresponding to the trace of the stability border 
on the face (ACD). For mi 3> = 3> m2, we have a heavy stable (1,3,4) ion which 
attracts the charge m-2, a stable situation according to the discussion in Sec. [VI]. 

The middle between A and C corresponds to masses m^ oc (1, 0, 1, 0), i.e. a configuration 
(pe + pe~), which seems unstable. We suspect there is no stability at all along AC, and 
similarly along AD, BC and BD. This guess has of course to be checked. 

On the other hand, the middle between A and B describes the hydrogen molecule 
(ppe~~e~) with m 8 = (1,1,0,0). This is a region of confortable stability. We guess that 
stability holds all the way along AB. There is no doubt (pp'e~e~) is stable as long as both 
m(p) ^> m(e) and m(p') ^> m(e), since the Born-Oppenheimer approximation holds in such 
case. Stability should remain, in our opinion, for m(p) 3> m(e) and m{p') > m(e), as long as 
the equality m 3 = m 4 is kept. This, indeed , suffices to make the two threshods (l,3) + (2,4) 
and (l,4) + (2,3) degenerate, and again, it is a general, though empirical, observation that 
competing thresholds often generate attraction; one could also say that the (pe~) and (p'e~) 
atoms have comparable Bohr radii, and their overlap might lead to some favourable exchange 
forces between them. 



The face (BCD) corresponds to mi = 0. As discussed in Sec. [\T|, the stability domain is 
generally dictated by the stability of the 3-body system (m^mg mj), which is described in 
Fig. Exceptions are found near the edges of (BCD), where we have more than one light 
particle. The situation is of course identical on the other faces. 

Consider now the tetrahedron (ABCD) with normalised inverse distances, in Fig. [|. The 
summit A corresponds to the {e + ppp), when one goes to A along the symmetry axis. This 
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is a stable protonium ion surronded by an electron, i.e. a stable 4-body configuration. The 
shape of the stability domain near A is more precisely seen on a cross-section (bed) parallel 
to (BCD) and close to A: the trace on (bed) looks exactly as the stability plot of Fig. [2|, 
because we are in the regime of of Sec. [VT| where one of the masses is much lighter than of 
the others. 

The sides AB and CD correspond to inverse masses (0, 0, x, 1 — x), with < x < 1, i.e. 
to configurations (ppee 1 ) with m(e) <C m(p) and m(e') -C rn(p), but the ratio of m(e) to 
m(e') arbitrary. We believe they are stable. 

The sides AC, AD, BC, and BD have inverse masses of the type (x, 0, 1 — x,0), i.e. 
(e + pe'p), for which dissociation seems likely to occur spontaneously. 

These speculations suggest the topology schematically drawn in Fig. [|. Note that with 
these variables, a face such as (BCD) describes situations of the type (pe + (e')~(e")~), where 
(with obvious notations) m, m', and m" are much lighter than M. It is stable for m = m' = 
m", with little excess binding We thus expect the stability band to be rather narrow 
near the centre of (BCD). 

VIII. LARGER MOLECULES 

An obvious extension of the present study consists of considering N = N + + iV_ > 4 
unit charges. The problem is in general rather complex, due in particular to the prolifer- 
ation of competing thresholds, but one can retain from the study in the N = 4 case that 
simplifications occur when very different masses are involved. 

Consider for instance N = 5, with masses and charges (M + M + M~m~m~), and M ^> m. 
The heavy core (M + M + M~) is stable, by simple rescaling of the positronium-ion case 
(e + e~e~). It acts as a localised positive charge, and should bind two electrons, with a 
wave function very similar to that of the familiar H~(pe~e~) ion. a similar reasoning was 
applied to molecular complexes involving muons and electrons JjlJ. If (M + M + M~m~m~) is 
stable for large M/m, a cooled antiproton could be trapped in a hydrogen molecule, before 
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annihilating. According to the current analysis ||20|| , when a p is captured in hydrogen, 
it quickly expels the electrons by Auger emission; a protonium atom (pp) is formed in a 
highly excited state, and it rapidly decays toward states of low orbital momemtum, where its 
annihilates. Recently some events with delayed annihilation have been reported, stimulating 
some theoretical studies some metastable states with an antiproton, a nucleus and 

some electrons might be formed in rare occasions. The above (pppe~e~) compound is another 
new possibility. 

Similarly, for N = 6, and iV + = AL = 3, the state (ppe + pe~e~) is probably stable, due 
to the combined stability of the (ppp) ion and of the positronium hybride (pe + e~e~). 

The symmetric case of N = 6, (m + m + m + m~m~m~), is much more delicate. Stability 
seems excluded for actual electrons, due to the Pauli principle, but one can address the 
question for pions, ir + and 7r~, which are bosons. The lowest threshold is (m + m + m~mr) + 
(m + m~), so one has to show that the lowest energy E decreases faster than — N, already 
for small N. For large N, a behaviour E cx —A 7 / 5 has been identified |23|. 



IX. APPLICATIONS TO HADRON SPECTROSCOPY 

A simple, but phenomenologically successful, description of the hadron spectrum is the 
non-relativistic quark model. Mesons are quark-antiquark bound states, and baryons consist 
of three quarks. An interesting property is flavour independence. At first approximation, i.e. 
without spin forces and relativistic corrections, the potential is the same whatever quarks, 
u, d, s, c, or b, are bound together. This is reminiscent from atomic physics, where the very 
same — r _1 potential acts in hydrogen, positronium, protonium, etc. 

In quark models with flavour independence, one expects some of patterns observed in 
atomic physics, those which are due to the universality of the potential, and independent of 
its particular Coulombic shape. For instance, the convexity property mentioned in Sec. |TJ 
for the 2-body energies is satisfied, and the inequality 

QQ + qq< 2Qq (25) 
13 



has indeed been noticed []10| , and observed in the experimental spectrum ||. 

We have seen that among the 4-body molecules, (ppe~e~) was the most stable. It is 
not surprising that explicit four-quark calculations P,p4"f have led to the prediction that the 
exotic configuration {QQqq) becomes stable if the mass ratio m(Q)/m(q) is large enough, 



whereas the equal-mass case QQQQ seems unbound. In this field, one call "exotic" a 
state whose quantum numbers, in particular flavour, cannot be matched by ordinary quark- 
antiquark or three-quark structures. Other multiquarks are predicted on the basis of spin- 
dependent forces, but this is out of the scope of this paper. It is hoped that double-charm, 
or charm-and-beauty spectroscopy will reveal some good surprises. 



X. CONCLUSIONS 

Let us summarise our results about 4 unit-charge systems, which bear a sharply decreas- 
ing degree of rigour. 

We first recalled the beautiful proof by Hylleraas and Ore of the stability of the positron- 
ium molecule [Q]. It can be supplemented by a rigourous study of the ordering of the various 
possible thresholds @||. 

Once the stability of (e + e + e~e~) is established, that of the hydrogen-like molecules 
(M + M + m~m~) follows rigourously HJ5J]. This is essentially a consequence of the variational 
principle. 

Now, one can prove a minimal extension of the stability domain toward other directions 
like (M + m + M~m~) or (mim^m^m^), if one takes for granted a non-rigourous input, 
namely the actual binding energy of (e + e + e~e~), which results from numerical variational 
calculations. 

The above results are far from covering all possible configurations. To picture how 
the stability domain could look like, we have substituted some missing results by simple 
guesses, which have to be checked. For instance, we feel that the asymmetric hydrogen 
molecule (pp'e~e~) remain stable whatever mass is given to p', as long as m(p) ^> m(e). 
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We have also assumed without proof that (pe + p'e~) immediately breaks into (pp') + (e + e~), 
except near the limit m{p') — > m(e). 

We also presented some plausible but non rigourous extensions to larger molecules and 
to multiquark systems in the quark model. 

In spite of its weaknesses, this study illustrates once more how quantum binding energies 
have a simple and well-controlled behaviour when one varies the constituent masses without 
changing the potential. It is hoped it will stimulate some interest for the investigations 
which remain to be done. 
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FIGURES 



FIG. 1. Graphical representation of the stability domain for 3 unit charges % = ±(— 1, 1, 1) 
with normalised masses m,, J2 m i = 1- The stability domain includes the symmetry axis (dotted 
line) where ni2 = m^. 



FIG. 2. Graphical representation of the stability domain for 3 unit charges qi = ±(—1, 1, 1) 
with normalised inverse masses ui = 1/mj, J2 a i = 1- The stability domain includes the symmetry 
axis (dotted line) where ct2 = 03. 
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FIG. 3. Schematic picture of the stability domain for 4 unit charges qi = (—1,-1,1,1). 
The stability frontier is shown by its trace on the unfolded faces of the tetrahedron (ABCD) of 
normalised masses rrii, J2 m i = 1- 
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FIG. 4. Schematic picture of the stability domain for 4 unit charges % = (—1, —1, 1, 1). The 
stability frontier is shown by its trace on the unfolded faces of the tetrahedron (ABCD) of inverse 
masses oti = 1/rrii, with normalisation J2 a i = 1- The trace of the stability domain in the section 
(bed) below the summit A looks like the 3-body domain of Fig. |2[ 
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